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ABSTRACT 


Boost  phase  interception  of  ballistic  missiles  is  envisioned  as  the  primary  re¬ 
sponse  of  the  layered  defense  architecture  implemented  in  the  ballistic  missile  defense 
system.  A  limited  time  frame  in  which  to  take  action  and  the  necessity  to  implement  hit- 
to-kill  technology  in  the  kill  vehicle  counterbalances  the  many  advantages  of  boost  phase 
interception.  Direct  hit  missile  technology  is  constrained  by  the  requirement  to  minimize 
miss  distance  to  a  negligible  amount  between  the  kill  vehicle  and  optimum  aimpoint  on 
the  target.  This  thesis  examines  kill  vehicle  effectiveness,  which  is  tantamount  to  miss 
distance,  as  a  function  of  both  the  kill  vehicle  maximum  acceleration  capability  and  the 
guidance  system  time  constant  necessary  to  destroy  a  target.  The  kill  vehicle  guidance 
system  is  modeled  in  MATLAB  as  a  fifth-order  binomial  series  with  proportional  naviga¬ 
tion.  The  simulation  examines  the  effect  of  an  accelerating  target  attributed  to  powered 
flight  and  aimpoint  displacement  caused  by  a  shift  in  tracking  point  from  the  target  plume 
to  the  payload  when  resolution  occurs.  The  kill  vehicle  minimum  requirements  as  indi¬ 
cated  by  the  simulation  include  a  lateral  acceleration  capability  of  four  times  the  target 
acceleration  and  a  guidance  system  time  constant  that  is  less  than  one -tenth  the  estimated 
flight  time. 
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EXECUTIVE  SUMMARY 


The  United  States  Department  of  Defense  is  developing  a  ballistic  missile  defense 
system  to  provide  a  layered  defense  capability  for  national  and  theater  level  protection. 
Boost  phase  interception  of  ballistic  missile  threats  is  envisioned  as  the  primary  response 
of  the  layered  defense  network.  Currently  operational  weapon  systems  can  engage  ballis¬ 
tic  missile  threats  in  the  later  stages  of  their  flight;  however,  there  is  a  lack  of  fielded  sys¬ 
tems  for  boost  phase  interception. 

A  limited  time  frame  in  which  to  take  action  and  the  necessity  to  implement  hit- 
to-kill  technology  in  the  kill  vehicle  counterbalances  the  many  advantages  of  boost  phase 
interception.  Direct  hit  missile  technology  is  constrained  by  the  requirement  to  minimize 
miss  distance  to  a  negligible  amount  between  the  kill  vehicle  and  optimum  aimpoint  on 
the  target  payload.  This  thesis  examines  kill  vehicle  effectiveness,  which  is  tantamount 
to  miss  distance,  as  a  function  of  both  the  kill  vehicle  maximum  acceleration  capability 
and  the  guidance  system  time  constant  necessary  to  destroy  a  given  target.  All  results 
promulgated  in  this  document  are  based  on  simulated  engagements  between  a  threat  bal¬ 
listic  missile  and  a  kill  vehicle  employing  proportional  navigation  guidance. 

The  kill  vehicle  minimum  requirements  include  a  lateral  acceleration  capability  of 
four  times  the  target  acceleration  and  a  guidance  system  time  constant  that  is  less  than 
one-tenth  the  estimated  flight  time  after  warhead  resolution.  Development  of  a  kill  vehi¬ 
cle  that  fulfills  these  requirements  is  one  step  towards  the  successful  implementation  of  a 
weapon  system  that  can  be  utilized  to  destroy  a  ballistic  missile  during  the  boost  phase  of 
its  flight. 
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I.  INTRODUCTION 


The  United  States  Department  of  Defense  is  developing  an  integrated  ballistic 
missile  defense  system  (BMDS)  to  provide  national  and  theater  level  protection  from  bal¬ 
listic  missile  attacks.  BMDS  will  incorporate  layered  defense  architecture  to  engage 
threat  missiles  at  the  earliest  opportunity.  Failure  to  eliminate  a  threat  will  initiate  addi¬ 
tional  engagements  by  different  weapon  systems  until  a  positive  kill  assessment  is  made. 
Currently,  there  are  fielded  weapon  systems  that  can  engage  a  ballistic  missile  during  the 
later  stages  of  its  flight;  however,  there  are  no  operational  systems  that  can  engage  a 
threat  missile  during  the  initial  stage  of  its  flight,  known  as  boost  phase. 

Various  weapon  systems  have  been  proposed  for  boost  phase  interception  (BPI) 
of  ballistic  missiles,  including  airborne  lasers,  space-based  intercept  missiles,  and 
ground-based  intercept  missiles.  These  systems  will  have  different  capabilities  and  limi¬ 
tations  and  shall  require  advanced  technologies  to  implement. 

This  thesis  is  applicable  to  ground-based  intercept  missiles  and,  in  particular,  ex¬ 
amines  the  requirements  for  kill  vehicle  effectiveness  during  boost  phase  interception  of 
ballistic  missiles.  The  kill  vehicle  is  the  payload  on  the  intercept  missile  and  may  be  con¬ 
sidered  an  armored  fist,  which  imparts  a  destructive  force  on  the  threat  missile.  Kill  ve¬ 
hicle  effectiveness  is  gauged  by  the  amount  of  miss  distance  at  termination  of  the  en¬ 
gagement.  Successful  boost  phase  interception  is  tantamount  to  complete  destruction  of 
the  target  and  necessitates  a  negligible  miss  distance  from  the  optimum  aimpoint  on  the 
target  warhead.  The  kill  vehicle  incorporates  a  guidance  system  to  reduce  the  miss  dis¬ 
tance  to  an  amount  suitable  for  complete  destruction  of  target.  During  the  boost  phase,  a 
target  ballistic  missile  is  characterized  by  a  monotonically  increasing  acceleration  rate 
attributed  to  powered  flight  and  a  brilliant  exhaust  plume.  The  effect  of  target  accelera¬ 
tion  on  miss  distance  is  analyzed  during  the  terminal  maneuvers  of  the  kill  vehicle,  to  de¬ 
termine  the  maximum  acceleration  requirement  for  a  kill  vehicle  to  impact  a  specific  tar¬ 
get.  Additionally,  the  effect  of  aimpoint  displacement  associated  with  transferring  the 
guidance  aimpoint  from  plume  to  warhead,  at  such  time  as  resolution  occurs,  is  exam¬ 
ined.  This  analysis  leads  to  determination  of  a  guidance  system  time  constant,  which  tol- 
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erates  an  aimpoint  shift  while  realizing  a  negligible  miss  distanee  at  termination  of  en¬ 
gagement.  The  analysis  is  based  on  a  MATLAB  simulation  of  the  kill  vehiele  guidanee 
system,  which  incorporates  a  fifth-order  binomial  series  representation  of  the  system  dy¬ 
namics  and  proportional  navigation  for  the  guidance  law. 

In  Chapter  II,  the  thesis  commences  with  a  review  of  layered  defense  and  its  ap¬ 
plicability  to  BMDS.  Integral  to  the  layered  defense  scheme  are  the  different  phases  of  a 
ballistic  missile’s  flight  trajectory;  consequently,  they  are  defined  in  this  chapter.  The 
final  section  identifies  the  advantages  and  disadvantages  of  boost  phase  interception. 

Chapter  III  provides  a  description  of  various  technologies,  which  must  be  devel¬ 
oped  or  further  refined  in  order  to  deploy  an  operational  kill  vehicle  for  boost  phase  in¬ 
terception  of  ballistic  missiles.  Direct  hit  weapons  technology  is  identified  as  the  neces¬ 
sary  means  of  conveying  a  destructive  force  in  the  BPI  scenario.  Furthermore,  sensors 
will  play  an  indispensable  role  in  the  successful  employment  of  any  weapon  system  in¬ 
tended  for  this  application.  A  possible  configuration  for  the  kill  vehicle  sensor  suite  is 
provided  in  this  section.  Lastly,  the  main  thrust  of  the  thesis  revolves  around  the  guid¬ 
ance  system  utilized  in  the  kill  vehicle;  therefore,  various  guidance  techniques  are  identi¬ 
fied,  and  it  is  reasoned  why  only  proportional  navigation  was  considered  suitable  for  this 
application.  A  review  of  ideal  proportional  navigation  is  presented  here  for  comparison 
purposes.  The  chapter  concludes  with  a  description  of  a  more  realistic  model  for  a  pro¬ 
portional  navigation  guidance  system,  which  utilizes  a  fifth-order  binomial  series  repre¬ 
sentation  for  the  system  dynamics  and  is  employed  in  the  ensuing  simulation. 

Chapter  IV  encompasses  the  simulation.  The  Runge-Kutta  numerical  technique  is 
identified  as  the  means  of  approximating  the  differential  equations  that  make-up  the 
model.  Subsequently,  the  system  dynamics  of  the  model  are  described.  The  binomial 
series  representation  of  the  system  dynamics  is  further  reduced  into  components  identi¬ 
fied  in  the  MATLAB  simulation.  The  results  are  sub-categorized  in  sections  for  accelera¬ 
tion  saturation  and  aimpoint  displacement. 

The  concluding  chapter  summarizes  the  results.  The  Appendices  include  a  listing 
of  the  MATLAB  code  used  in  the  simulation,  and  a  derivation  of  the  seeker  transfer  func¬ 
tion  used  in  the  guidance  simulation. 
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II.  BALLISTIC  MISSILE  DEFENSE 


Prior  to  January  2002,  the  United  States  missile  defense  program  was  organized 
within  two  frameworks,  national  missile  defense  (NMD)  and  theater  missile  defense. 

The  former  was  intended  to  protect  the  continental  United  States  from  attacks  by  inter¬ 
continental  ballistic  missiles  (ICBM)  while  the  latter  was  intended  to  protect  US  military 
forces  deployed  overseas  against  all  missile  attacks,  up  to  and  including  intermediate 
range  ballistic  missiles  (IRBM). 

In  January  2002,  the  Secretary  of  Defense  ordered  the  establishment  of  a  single 
integrated  program  to  develop  a  ballistic  missile  defense  system  capable  of  defending  the 
United  States,  its  deployed  forces,  and  allies  [1].  BMDS  shall  implement  a  layered  de¬ 
fense  battle  doctrine,  capable  of  intercepting  ballistic  missiles  in  all  phases  of  their  flight 
trajectory.  This  chapter  will  describe  the  principle  of  layered  defense  and  its  application 
to  ballistic  missile  defense. 

A,  LAYERED  DEFENSE 

The  principle  of  layered  defense  has  various  interpretations;  it  is  instructive  to 
first  examine  the  naval  combat  version,  which  is  a  precursor  to  the  model  utilized  by  the 
Missile  Defense  Agency  (MDA).  In  naval  operations,  threat  evaluation  and  weapons  as¬ 
signment  is  based  on  range,  bearing,  and  range  rate  of  the  target  to  the  defending  unit.  If 
a  positive  kill  assessment  is  not  made  after  the  primary  engagement,  the  defending  unit 
progressively  assigns  different  weapons  systems  for  target  engagement,  based  on  the 
range  from  the  defending  asset  to  the  target.  The  Missile  Defense  Agency  has  juxtaposed 
this  concept  with  the  different  segments  of  a  ballistic  missile’s  flight  trajectory.  The  lay¬ 
ered  defense  scheme  put  forth  by  MDA  [2]  includes  the  following  segments  in  sequential 
order: 

•  Boost  Segment; 

•  Midcourse  Segment;  and 

•  Terminal  Segment. 

Each  phase  of  the  missile’s  trajectory  is  associated  with  a  specific  altitude  sector 
above  the  earth’s  surface,  as  depicted  in  Figure  1.  MDA  intends  to  engage  and  destroy 
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the  threat  missile  at  the  earliest  possible  opportunity;  however,  if  this  is  unsuccessful,  a 
handover  will  be  made  to  the  next  weapon  system  in  the  layered  defense  chain. 


LEGEND 

LOWER  TIER 
ENDO-ATMOSPHERIC 
BOOST  PHASE  & 

TERMINAL  PHASE  (NOT  SHOWN) 

BOUNDARY  LAYER 
ENDO/EXO-ATMOSPHERIC 
MIDCOURSE  PHASE 

UPPER  TIER 
EXO-ATMOSPHERIC 
MIDCOURSE  PHASE 


DISTANCE  FROM  DEFENDING  ASSET 


Figure  1.  Layered  Defense  (After  Ref.  2). 


B.  BALLISTIC  MISSILE  TRAJECTORY 

A  ballistic  missile  trajectory  is  comprised  of  the  following  phases:  boost  phase, 
midcourse  phase,  and  terminal  phase  as  shown  in  Figure  2.  Each  of  these  phases  is  fur¬ 
ther  discussed  in  the  following  subsections. 


Figure  2.  Ballistic  Missile  Trajectory  (From  Ref.  3). 
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1,  Boost  Phase 

The  boost  phase  is  initiated  with  roeket  ignition  just  prior  to  liftoff  and  completed 
after  burnout  of  the  final  stage  of  the  ballistic  missile.  During  this  time  frame,  the  ballis¬ 
tic  missile  accelerates  to  its  maximum  speed.  The  missile  velocity  is  relatively  slow  dur¬ 
ing  this  phase  because  the  initial  velocity  is  zero  and  the  rocket  engine’s  propulsive  force 
is  opposed  by  the  earth’s  gravity.  The  discharge  of  hot  propellant  gasses  makes  the  ex¬ 
haust  plume  highly  visible  to  infrared  detectors  at  great  ranges.  Thrust  is  terminated  at 
the  end  of  boost  phase. 

The  duration  of  boost  phase  is  dependent  on  the  propellant  burn  rate,  the  missile 
range  capability,  and  the  average  acceleration  as  shown  in  Figure  3.  The  line  graph  at¬ 
tributed  to  an  acceleration  of  3g  is  indicative  of  older  missile  types,  such  as  the  SCUD, 
No  Dong  1  (ND-1),  and  the  Chinese  Surface-to-Surface  Missile  2  (CSS-2);  whereas 
newer  missiles  have  an  average  acceleration  of  approximately  lOg. 


Range  (km) 


Figure  3.  Boost  Phase  Bum  Times  for  IRBM  (After  Ref.  2). 
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Intercontinental  ballistic  missiles  have  a  range  capability  in  excess  of  5000  kilo¬ 
meters  and  boost  phase  duration  of  approximately  240  seconds  for  rocket  engines  using 
liquid  propellant  and  180  seconds  for  solid  propellant  rocket  engines  [4], 

ICBMs  will  exit  the  earth’s  atmosphere  by  the  end  of  boost  phase  whereas  inter¬ 
mediate  range  ballistic  missiles  only  reach  the  periphery  of  the  exo-atmosphere.  The 
maximum  altitude  and  downrange  attained  by  the  end  of  boost  phase  is  dependent  on 
various  design  parameters  including  the  guidance  law  utilized,  propellant  bum  rate,  ac¬ 
celeration,  and  launch  angle.  Typical  values  are  listed  in  Table  1. 


Missile  Range  (km) 

Burnout  Altitude  (km) 

Burnout  Range  (km) 

500 

20-40 

25-75 

1,000 

45-70 

40-100 

2,000 

70-130 

75-150 

3,000 

100-170 

125-250 

10,000 

175-220 

425-475 

Table  1.  Approximate  Coordinates  at  Termination  of  Boost  Phase  (After  Ref  2). 

2.  Midcourse  Phase 

The  midcourse  phase  commences  after  termination  of  the  boost  phase  and  en¬ 
compasses  the  longest  segment  of  the  trajectory,  approximately  20  to  30  minutes  for  an 
ICBM  [3].  The  weapon  payload,  which  may  be  a  unitary  warhead,  multiple  re-entry  ve¬ 
hicles,  or  submunitions,  is  deployed  during  this  phase  and  follows  an  unpowered  ballistic 
trajectory  towards  the  intended  target.  The  midcourse  trajectory  will  take  place  in  the 
exo-atmosphere  for  ICBMs  and  intermediate  range  ballistic  missiles. 

Intercepting  the  weapon  payload  during  the  midcourse  phase  is  complicated  by 
the  need  to  eliminate  all  submunitions  or  re-entry  vehicles  (RV),  plus  discriminating  and 
excluding  any  countermeasures  and  booster  debris.  Figure  4  provides  a  pictorial  repre¬ 
sentation  of  the  crowded  scenario  encountered  during  the  midcourse  phase;  also  shown  is 
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the  attitude  control  module  (ACM),  which  corrects  small  guidance  errors  and  ensures  the 
weapon  maintains  the  desired  course. 


Figure  4.  Unopposed  Ballistic  Missile  Flight  (From  Ref.  5). 


3,  Terminal  Phase 

The  terminal  phase  is  the  ending  segment  of  the  missile’s  flight  trajectory.  The 
weapon  payload  has  amassed  great  speed  by  the  start  of  this  phase;  however,  re-entry  into 
the  earth’s  atmosphere  will  impart  changes  to  the  weapon’s  velocity  and  trajectory.  At¬ 
mospheric  drag  will  ensure  the  weapon  is  separated  from  the  debris  field.  Lastly,  the 
weapon  is  capable  of  maneuvering  to  avoid  interception  during  this  phase. 

C.  BOOST  PHASE  INTERCEPTION 

Boost  phase  interception  of  a  ballistic  missile  threat  shall  be  the  primary  response 
of  the  BMDS  layered  defense  architecture.  Boost  phase  interception  has  many  advan¬ 
tages  and  disadvantages,  identified  in  the  following  subsections  [6]. 

1,  Advantages  of  Boost  Phase  Interception 

The  advantages  of  a  response  in  this  segment  of  the  layered  defense  scenario  in¬ 
clude  the  following; 

•  Single  target  engagement  preempts  deployment  of  multiple  re-entry  vehi¬ 
cles  or  submunitions  from  threat  missile  payload; 

•  Large  infrared  signature  attributed  to  the  exhaust  plume  makes  ballistic 
missile  readily  detectable  at  long  ranges; 

•  Low  velocity  of  threat  missile  at  commencement  of  boost  phase  permits 
employment  of  terrestrial  interceptor; 
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•  Forward  defense  presents  opportunity  to  destroy  ballistic  missile  far  away 
from  the  United  States; 

•  Sensors  utilized  for  BPI  will  serve  as  an  early  warning  mechanism  for  the 
additional  defensive  layers;  and 

•  Any  ballistic  missiles  that  penetrate  this  outer  layer  will  be  identified  for 
handover  and  weapons  assignment  to  the  inner  layers  of  the  defense  net¬ 
work. 

2,  Disadvantages  of  Boost  Phase  Interception 

Boost  phase  interception  is  constrained  by  the  following: 

•  Limited  time  span  in  which  to  engage  target  necessitates  development  of 
advanced  weapons  systems  to  accomplish  this  mission; 

•  A  ground-based  intercept  system  must  be  located  relatively  close  to  launch 
point  of  threat  missile; 

•  The  variable  acceleration  rate  attributed  to  fuel  consumption  and  discard¬ 
ing  of  spent  stages  complicates  targeting; 

•  BMDS  sensors  cannot  pre-determine  total  pitch/rotation  the  missile  must 
undergo  to  attain  ballistic  trajectory; 

•  Advanced  ballistic  missiles  can  maneuver  for  energy  management  or  in¬ 
tercept  avoidance; 

•  Debris  from  the  engagement,  including  undestroyed  submunitions,  may 
fall  on  friendly  territory;  and 

•  The  kill  vehicle  (KV)  must  incorporate  advanced  hit-to-kill  technology  to 
ensure  complete  destruction  of  the  target,  including  all  submunitions. 


The  development  of  an  effective  kill  vehicle  for  boost  phase  interception  is  a  ma¬ 
jor  undertaking,  which  necessitates  the  integration  of  sensors,  lethality  enhancement 
mechanisms,  and  guidance  and  control  systems  in  a  payload  that  is  deployable  from  an 
interceptor  missile.  A  description  of  these  technologies  is  provided  next  in  Chapter  III. 
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III.  KILL  VEHICLE  TECHNOLOGIES 


During  the  Cold  War  era,  missile  interceptors  armed  with  conventional  explosive 
warheads  were  deemed  ineffective  at  countering  ballistic  missile  threats.  Consequently, 
the  United  States  and  the  former  Soviet  Union  developed  nuclear  armed  interceptors  for 
employment  as  anti-ballistic  missiles  [7].  This  solution  was  politically  unacceptable  since 
a  nuclear  explosion  at  high  altitudes  over  a  defended  country  would  incapacitate  numer¬ 
ous  electrical  and  electronic  systems  in  that  country.  The  current  solution  is  implementa¬ 
tion  of  direct  hit  or  hit-to-kill  technology,  whereby  the  impact  of  collision  is  sufficient  to 
destroy  the  target.  The  primary  impediment  of  this  technology  is  the  requirement  for  in¬ 
finitesimal  miss  distance  between  the  kill  vehicle  and  optimum  aimpoint  on  the  target 
warhead,  to  ensure  destruction  of  all  submunitions  [8].  This  chapter  considers  direct  hit 
weapons  technology,  the  sensor  suite,  and  guidance  system  requirements  for  a  kill  vehicle 
suitable  for  boost  phase  interception. 

A,  DIRECT  HIT  WEAPON  TECHNOLOGY 

The  destructive  force  conveyed  by  a  direct  hit  weapon  is  several  orders  of  magni¬ 
tude  greater  than  a  conventional  blast  fragmentation  warhead. 

This  focused  direct  hit  impact  energy  is  created  by  a  nearly  solid  concen¬ 
tration  of  directed  energy,  but  must  rely  on  the  radial  expansion  of  that  en¬ 
ergy  to  damage  any  submunitions  not  located  on  the  impacting  missile 
flight  path  [8]. 

A  major  constraint  of  direct  hit  weapons  is  locating  the  aimpoint  on  the  target  payload 
where  the  maximum  number  of  submunitions  will  be  destroyed.  This  optimum  aimpoint 
is  euphemistically  known  as  the  sweet  spot.  The  lethality  of  a  direct  hit  weapon  will  de¬ 
teriorate  significantly  as  the  aimpoint  is  displaced  from  the  sweet  spot.  This  is  attributed 
to  decreased  radial  energy  imparted  to  the  area  offset  from  the  aimpoint. 

The  ALPHA-KV  analytical  model  developed  by  Orphal  [9]  graphically  illustrates 
that  aimpoint  errors  necessitate  increased  kinetic  energy  for  the  kill  vehicle  to  achieve 
lethality.  The  model  assumes  that  the  crater  volume  ,  caused  by  the  kill  vehicle  colli¬ 
sion  with  the  target,  is  proportional  to  the  impact  kinetic  energy 
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where  m  is  the  kill  vehiele  mass  and  V  is  the  elosing  veloeity.  The  proportionality  eon- 
stant  a  (em  /J)  is  determined  experimentally  or  by  using  a  cratering  model  such  as  the 
Tate  model  [10].  The  Tate  theory  of  penetration  states  that  a  can  be  approximated  by 
the  relationship 
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where  Pj  and  are  the  bulk  densities  of  the  target  and  kill  vehicle,  respectively.  The 
parameter  Rj  (N/cm  )  is  determined  from  the  sound  speed  of  the  target  and  the 
Brinell  hardness  number  (BHN)  of  the  target  in 


=4.2x1  O' BHN 


3  ^  4.2x10'  ) 


The  required  crater  volume  is  assumed  cylindrical  and  is  given  by 

Vy,>nLc[Ro^R„) 


(3) 


(4) 


where  is  the  crater  depth,  R^  is  the  minimum  crater  radius  to  overlap  all  submunitions, 
and  is  any  offset  from  the  optimum,  as  illustrated  in  Figure  5  where  the  submunitions 
are  highlighted  in  gray.  For  maximum  effectiveness,  which  is  tantamount  to  zero  miss 
distance,  the  optimum  hit-point  would  be  the  geometric  center  of  the  payload.  If  the  hit- 
point  is  offset  by  an  amount  R^^  from  the  optimum  location  the  crater  radius  will  have  to 

be  increased  to  an  amount  equivalent  to  R^  +  K- 
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Figure  5.  Target  Payload  with  Optimum  Hit-Point  (After  Ref.  9). 

Ineorporating  target  parameters  in  the  model  ealculations,  ineluding  number  of 
submunition  tiers,  diameter  of  individual  submunitions,  and  overall  length  and  diameter 
of  target,  optimizes  the  ALPHA-KV  model.  The  model  is  easily  programmed  in 
MATLAB,  attaehed  as  Appendix  A,  and  is  used  to  generate  Figure  6  whieh  illustrates 

that  the  kinetie  energy  [mV^ j 2^  required  by  the  kill  vehiele  is  near  linearly  proportional 
to  the  amount  of  offset  )  from  the  sweet  spot. 

The  exaet  loeation  of  the  sweet  spot  is  unique  to  eaeh  missile  type  depending  on 
the  payload  and  geometry  and  eannot  be  aseertained  preeisely  with  eurrent  sensing  sys¬ 
tems;  therefore,  a  default  loeation  must  be  utilized.  This  default  loeation  must  be  suitable 
for  all  threats,  whieh  the  direet  hit  weapon  is  designed  to  engage  and  destroy.  Loeation 
of  the  sweet  spot  is  not  as  eritieal  with  respeet  to  targeting  unitary  warheads. 
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Figure  6.  Required  Kinetic  Energy  for  Optimum  Kill  Vehicle. 


Kill  vehicles  for  boost  phase  interception  have  not  been  fully  implemented  to 
date;  however,  similar  systems  for  ground-based  midcourse  defense  have  been  designed 
and  tested.  Raytheon  has  developed  and  successfully  trialed  the  exo-atmospheric  kill  ve¬ 
hicle  (EKV)  pictured  in  Eigure  7.  Visible  in  this  schematic  representation  are  two  of  four 
thrusters  arranged  in  a  cruciform  configuration  and  used  to  maneuver  the  KV  on  the  hori¬ 
zontal  (pitch)  axis  and  vertical  (yaw)  axis.  The  thrusters  are  part  of  the  divert  and  attitude 
control  system  (DACS)  which  steers  the  kill  vehicle  in  accordance  with  the  guidance 
commands.  Most  of  the  forward  momentum  required  by  this  kill  vehicle  is  provided  by 
the  intercept  booster’s  burnout  velocity.  The  kill  vehicle  weighs  approximately  140 
pounds  and  has  a  length  of  55  inches  and  a  diameter  of  24  inches  [3].  Also  identifiable  in 
the  left  foreground  of  Eigure  7  is  the  acquisition  telescope.  The  sensor  system  incorpo¬ 
rates  three  separate  sensors,  not  shown  in  Eigure  7,  which  measure  visible  light,  mid¬ 
wave  infrared,  and  long-wave  infrared  radiation. 
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Figure  7.  Exo-atmospheric  Kill  Vehicle  (From  Ref.  5). 


B.  SENSOR  SUITE 

The  principal  signature  of  a  ballistic  missile  during  its  boost  phase  trajectory  is  at¬ 
tributed  to  the  exhaust  plume.  It  is  characterized  by  strong  emissions  in  the  mid-wave 
infrared  (MWIR)  spectrum,  3  to  5  jum ,  and  weaker  emissions  in  the  ultraviolet  spectrum 
[11].  Another  valuable  source  of  radiation  is  skin  friction  from  the  missile  hardbody, 
which  radiates  strongly  in  the  long-wave  infrared  (LWIR)  band,  8  to  12  jum . 

Launch  detection  and  tracking  of  a  ballistic  missile  threat  will  be  conducted  by  a 
space  based  infrared  (IR)  system,  which  will  sense  the  target  exhaust  plume  and  provide 
a  two-dimensional  angular  measurement  of  bearing  and  elevation.  Ground-based  air  de¬ 
fense  radars  will  acquire  and  track  the  target  after  it  climbs  above  the  radar  horizon. 

They  will  provide  a  three  dimensional  measurement  of  target  range,  bearing,  and  eleva¬ 
tion.  The  two  data  types,  infrared  and  radio  frequency,  will  be  fused  into  one  threat  track 

'i 

and  passed  to  the  Battle  Management/Command,  Control,  and  Communications  (BM/C  ) 
system  for  engagement  authority  [12].  If  permission  to  engage  is  granted,  the  target  will 
be  designated  for  weapons  assignment.  A  ground-based  intercept  (GBI)  missile  will  con¬ 
sequently  be  launched  at  the  threat  missile.  Various  proposals  for  intercept  launchers 
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have  been  made;  however,  this  study  will  only  examine  ground-based  interceptors.  It  is 
assumed  that  the  interceptor  will  incorporate  a  global  positioning  system  (GPS)  aided  in- 
ertial  navigation  system,  and  communication  link  with  the  BM/C  system.  Initially,  the 
interceptor  will  guide  to  the  planned  intercept  point  (PIP)  generated  from  the  fused  data. 
It  will  receive  target  updates  on  a  continuous  basis  until  the  kill  vehicle  is  deployed,  and 
its  onboard  sensors  takeover  the  sensing  mission. 

The  sensor  suite  onboard  the  kill  vehicle  should  include  a  dual-band  imaging 
seeker,  operating  in  the  MWIR  and  LWIR  bands.  It  is  implemented  with  two  co¬ 
registered  quantum- well  infrared  photodetector  (QWIP)  focal  plane  arrays  (FPA)  [13]. 
The  advantage  of  this  dual-band  IR  camera  is  simultaneous  imaging  of  the  rocket  plume 
and  missile  hardbody.  Digital  image  processing  (DIP)  is  required  for  comparison  of  the 
simultaneous  MWIR  and  LWIR  images.  The  DIP  algorithm  is  known  as  MWIR/LWIR 
band-ratioing.  A  plot  of  the  in-band  radiance  arriving  at  the  MWIR  detector  versus  that 
of  the  LWIR  detector  will  generate  a  line,  with  a  slope  equivalent  to  the  ratio  between  the 
two  bands.  Planck’s  blackbody  radiance  function  is  first  evaluated  as  a  precursor  to  de¬ 
termining  the  band  ratio.  The  radiance  function  is  given  by  [14] 

2hr^ 

A  [exp(Ac//lA:r)-lJ 


where  h  is  Planck’s  constant,  c  is  the  speed  of  light,  k  is  Boltzmann’s  constant,  A  is 
wavelength,  and  T  is  temperature.  The  in-band  radiance  incident  on  the  sensor  optics  is 
then  determined  by  integration  of  Planck’s  blackbody  radiance  function  with  respect  to 
the  wavelengths  of  interest.  The  ratio  can  be  evaluated  by 


j  B,{T)dA 

Ratio  =  ^ - . 

1  BAT)dl 

LWIR 


(6) 


A  proof  of  concept  for  dual-band  ratioing  was  developed  in  MATLAB  and  is 
available  from  the  thesis  advisor.  The  model  evaluates  Equations  (5)  and  (6)  and  is  used 
to  plot  the  numerator  versus  denominator  of  the  Ratio  as  shown  in  Figure  8. 


14 


2.2 


Figures.  MWIR/LWIR  Band-Ratioing. 

Figure  8  shows  unique  lines  are  generated  for  the  various  target-combinations  that  may 
be  in  the  sensor’s  field  of  view,  such  as  only  the  plume  at  great  distances,  or  both  the 
plume  and  missile  hardbody  as  the  KV  closes  the  target,  or  only  the  hardbody  during  the 
terminal  maneuvers  of  the  engagement. 

In  conjunction  with  band-ratioing,  background  subtraction  is  conducted  to  adjust 
the  gain  on  each  of  the  two  image  sensors,  MWIR  and  LWIR.  This  will  eliminate  back¬ 
ground  sources  of  radiation  such  as  the  sky  and  sun  to  ensure  that  the  band-ratioing 
calculations  are  only  evaluated  for  target  sources. 

Initially,  the  target  plume  will  produce  a  dominant  signature  in  the  MWIR  band. 
The  band-ratioing  technique  will  generate  values  similar  to  those  for  the  plume  only  line 
in  Figure  8  and  the  kill  vehicle  will  track  on  the  power  centroid  of  the  plume.  As  the  kill 
vehicle  closes  on  the  target,  skin  friction  emissions  associated  with  the  BM  hardbody  will 
become  highly  visible  in  the  LWIR  band.  The  band-ratioing  technique  will  generate  val¬ 
ues  similar  to  the  hardbody  only  line  in  Figure  8,  and  the  kill  vehicle  will  switch  to  track- 
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ing  the  leading  edge  of  the  missile  hardbody.  Resolution  of  the  hardbody  will  permit  dis¬ 
crimination  of  the  optimum  aimpoint  on  the  nosecone  section  of  the  ballistic  missile. 

This  is  accomplished  by  displacing  the  aimpoint  from  the  leading  edge  to  the  specified 
default  location. 

The  applicability  of  dual  band  imaging  is  demonstrated  in  Figure  9,  which  exhib¬ 
its  an  Atlas  V  rocket,  2  minutes  after  liftoff,  at  a  downrange  of  35  km  and  an  altitude  of 
24  km.  The  plume  signature  is  saturated  in  both  bands;  however,  the  missile  hardbody  is 
clearly  visible  in  the  LWIR  band. 


Figure  9.  Dual  Band  IR  Image  of  Atlas  V  Launch  21  August  2002  (From  Ref  13). 

C.  GUIDANCE 

Various  techniques  have  been  proposed  for  strategic  missile  guidance  including 
proportional  navigation  (PN),  augmented  proportional  navigation,  predictive  guidance, 
and  optimal  control  guidance  (OCG)  [6].  Proportional  navigation  is  prevalent  in  current 
missile  applications  and  considered  the  most  effective  guidance  technique  when  accurate 
scientific  or  technical  intelligence  is  lacking  about  the  target  missile.  For  example,  pre¬ 
dictive  guidance  requires  knowledge  of  the  threat  missile’s  dynamics  for  best  perform¬ 
ance,  and  OCG  systems  are  difficult  to  mechanize  and  have  not  been  implemented  to  date 
[15].  The  proportional  navigation  guidance  law  provides  a  basis  for  the  simulation  de¬ 
veloped  in  this  thesis. 

1.  Ideal  Proportional  Navigation 

The  definition  of  proportional  navigation  is  stated  as: 

The  rate  of  change  of  missile  heading  is  directly  proportional  to  the  rate  of 
rotation  of  the  line-of-sight  (LOS)  from  the  missile  to  the  target  [16]. 
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The  proportional  navigation  guidance  law  is  mathematically  expressed  as 

He  =  N'Vc  ^  (7) 

where  is  the  commanded  normal  acceleration,  N'  is  the  effective  navigation  ratio 
which  determines  the  extent  of  damping  in  the  associated  control  system,  is  the  clos¬ 
ing  velocity  between  kill  vehicle  and  target,  and  X  is  the  angular  line-of-sight  rate  from 
kill  vehicle  to  target  ballistic  missile.  The  effective  navigation  ratio  is  a  dimensionless 
parameter,  empirically  determined  to  be  most  effective  in  the  range  of  3  to  5 .  An  ideal 
representation  of  proportional  navigation  excludes  system  dynamics  and  is  known  as  a 
zero-lag  guidance  system  as  shown  in  Figure  10. 
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Figure  10.  Ideal  Representation  of  Proportional  Navigation  (After  Ref.  17). 


The  variables  displayed  in  Figure  8  include  the  achieved  lateral  acceleration  of  the 
kill  vehicle  ,  the  target  normal  acceleration  rij ,  the  commanded  acceleration  ,  and 


the  relative  acceleration  between  target  and  kill  vehicle  y  as  given  by 


y  =  n^-n^. 


(8) 
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The  relative  aeeeleration  is  transformed  by  double  integration  into  relative  position  y  , 
whieh  is  equal  to  the  miss  distance  at  the  end-of-flight  time  tp .  The  time-to-go  until  in¬ 
tercept  completion  is  given  by 

(9) 

where  t  is  the  current  time.  Division  of  the  relative  position  by  the  range-to-target 
{Vctgo)  produces  the  geometric  line-of-sight  angle /I  as  follows 
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The  seeker  tracks  the  target  by  differentiation  of  A.  to  generate  A  .  Idealized  system  dy¬ 
namics  permit  unity  gain  representation  of  the  noise  filter.  The  commanded  acceleration 
Hp,  as  previously  stated  in  Equation  1  is  indicative  of  proportional  navigation  and  pro¬ 
vides  the  guidance  commands  to  the  flight  control  system.  Lastly,  the  flight  control  sys¬ 
tem  employs  small  rocket  thrusters  to  maneuver  the  KV  with  respect  to  vertical  and  hori¬ 
zontal  axes  to  ensure  a  direct  hit  with  the  optimum  aimpoint  on  the  target.  The  flight  con¬ 
trol  system  is  also  designated  as  unity  gain  because  of  the  idealized  representation.  The 
achieved  lateral  acceleration  is  the  end  result  of  the  flight  control  system  responding 
to  the  guidance  commands. 

If  this  ideal  system  were  physically  realizable,  every  engagement  scenario  would 
result  in  a  constant-bearing  course  with  parallel  consecutive  line-of-sight  paths  and  as¬ 
sured  destruction  of  the  target,  assuming  acceleration  saturation  was  preempted.  In  other 
words,  the  miss  distance  would  always  be  reduced  to  zero  at  termination  of  the  engage¬ 
ment.  Consequently,  a  higher-order  model  of  the  guidance  system  dynamics  is  required 
to  preclude  an  unrealistic  representation  of  proportional  navigation  systems. 

2.  Binomial  Series  Representation  of  PN  Guidance  System  Dynamics 

A  binomial  series  representation  of  the  guidance  system  dynamics  provides  an  ef¬ 
fective  means  for  modeling  a  more  realistic  proportional  navigation  system  [17].  A  rep¬ 
resentation  of  the  system  dynamics  in  the  s  domain  is  given  by 
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where  r  is  the  guidance  system  time  constant,  n  is  the  system  order,  and  H.O.T.  repre¬ 
sents  higher-order  terms,  which  have  a  minor  effect  on  the  model  and  are  omitted  for 
simplicity.  A  fifth-order  representation,  considered  suitable  for  a  preliminary  analysis  of 
kill  vehicle  performance,  allocates  five  equivalent  time  constants  as  follows:  one  for  the 
seeker,  one  for  the  noise  filter,  and  three  for  the  flight  control  system  in  Figure  11. 


FLIGHT  CONTROL 
SYSTEM 


GUIDANCE 


NOISE  FILTER 


Figure  11.  Fifth-Order  Representation  of  PN  Guidance  System  (After  Ref.  18). 

In  this  representation,  the  noisy  line-of-sight  rate  must  be  filtered  to  produce  a 

better  estimate  of  X  .  Subsequently,  the  guidance  commands  are  generated  on  the  basis 
of  the  filtered  output.  In  this  model  of  PN,  the  response  to  the  guidance  commands  is 
also  affected  by  the  flight  control  system  dynamics.  Also  shown  in  Figure  1 1  is  the  pa¬ 
rameter  ,  which  is  used  to  enter  the  aimpoint  displacement  as  an  initial  condition  of 

the  guidance  simulation.  The  transfer  function  for  the  system  shown  in  Figure  9  is  given 
by 
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The  effect  of  system  dynamics  on  the  calculation  of  miss  distance  is  made  evident 
in  Figure  12,  which  shows 

how  the  miss  distance  due  to  a  3g  target  maneuver  varies  with  flight  time 
and  system  order  for  a  binomial  guidance  system  in  which  the  effective 
navigation  ratio  is  4  and  the  effective  guidance  system  time  constant  is  1  s. 

We  can  see  that  the  performance  projections  resulting  from  a  single-lag 
guidance  system  model  are  a  serious  underestimate  of  the  influence  of  tar¬ 
get  maneuver  on  miss  when  the  flight  time  is  not  an  order  of  magnitude 
greater  than  the  guidance  system  time  constant  [17]. 


Figure  12.  Comparison  of  System  Dynamics  (From  Ref.  17). 


It  is  apparent  in  Figure  12  that  the  first-order  system  will  underestimate  the 
amount  of  miss  distance  generated.  For  example,  in  an  engagement  of  duration  3  sec¬ 
onds,  the  first-order  system  indicates  a  negligible  miss  distance,  whereas  the  fifth-order 
system  indicates  a  miss  distance  of  45  feet.  If  these  engagements  transpired  in  actuality, 
the  kill  vehicle  with  a  first-order  system  would  assume,  just  prior  to  termination  of 
engagement,  that  it  was  going  to  hit  the  aimpoint  and  cease  to  make  corrections.  The  kill 


20 


vehicle  with  the  fifth-order  system  would  continue  to  make  corrections  to  reduce  the  miss 
distance.  The  negative  segment  of  each  curve  indicates  that  the  kill  vehicle  is  below  the 
target  while  the  positive  segment  is  above  the  target.  Based  on  the  greater  accuracy  pro¬ 
vided  by  the  fifth-order  system,  it  was  used  in  the  ensuing  simulation. 

The  operability  of  a  kill  vehicle  for  BPI  is  predicated  on  the  integration  of  the 
various  technologies  described  in  this  chapter.  To  accomplish  the  seemingly  impossible 
task  of  hitting  a  bullet  with  a  bullet,  as  some  pundits  have  described  this  undertaking,  the 
efficacy  of  the  guidance  and  control  system  is  paramount.  The  next  chapter  presents  a 
simulation  of  the  kill  vehicle  guidance  system  operating  in  a  boost  phase  intercept  envi¬ 
ronment.  The  simulation  enables  determination  of  the  acceleration  required  by  the  kill 
vehicle  and  a  suitable  time  constant  for  the  guidance  system  dynamics. 
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IV.  SIMULATION 


The  simulation  models  the  terminal  maneuvers  of  the  kill  vehiele,  following  sepa¬ 
ration  from  the  last  stage  of  the  intereept  booster.  Prior  to  eommeneement  of  the  simula¬ 
tion,  it  is  assumed  that  the  kill  vehiele  has  target  loek-on  and  is  traeking  the  ballistie  mis¬ 
sile  exhaust  plume  via  radiant  emissions  in  the  MWIR  band.  Until  sueh  time,  the  kill  ve¬ 
hiele  guides  to  an  aimpoint  at  the  power  eentroid  of  the  target  exhaust  plume.  As  the  kill 
vehiele  eloses  on  the  target,  the  ballistie  missile  hardbody  will  beeome  visible  in  the 
LWIR  band.  At  this  time,  the  ballistie  missile  warhead  is  resolved,  and  the  guidanee  aim 
point  shifts  almost  instantaneously  from  the  plume  to  the  warhead.  This  apparent  target 
displaeement  oeeurs  late  in  the  intereept  seenario  and  is  modeled  as  a  step  input  to  the  kill 
vehiele  guidance  system.  The  simulation  demonstrates  the  affect  of  aimpoint  displace¬ 
ment  on  miss  distance  and  is  used  to  determine  a  time  constant  suitable  for  the  guidance 
system  dynamics  to  reduce  miss  distance  to  a  negligible  amount. 

It  has  been  previously  noted  that  a  ballistic  missile  accelerates  for  the  duration  of 
boost  phase;  therefore,  the  simulation  is  also  used  to  determine  what  acceleration  capabil¬ 
ity  is  required  by  the  kill  vehicle  to  preclude  saturation  and  ensure  lethal  impact  with  the 
accelerating  target. 

The  original  simulation  was  developed  by  Zarchan  [17]  to  study  the  effects  of  an 
apparent  target  displacement  that  occurs  when  the  sensor  tracking  point  shifts  from  the 
power  centroid  of  multiple  unresolved  targets  to  one  resolved  target. 

A,  NUMERICAL  METHOD 

The  kill  vehicle  guidance  system  is  modeled  as  a  fifth-order  binomial  series  repre¬ 
sentation  in  MATLAB.  The  model  utilizes  differential  equations  to  represent  the  in¬ 
put/output  relationships  of  the  seeker,  noise  filter,  and  flight  control  system. 

A  second-order  Runge-Kutta  numerical  integration  method  is  used  to  solve  the 
differential  equations  describing  the  kill  vehicle  guidance  system  [19].  The  main  diffi¬ 
culty  in  the  numerical  solution  of  ordinary  differential  equations  is  the  solution  of  the 
first-order  equation; 
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(13) 


dt 


with  respect  to  the  initial  condition  of  y  =  at  t  =  t^.  The  Runge-Kutta  technique  gen¬ 
erates  approximations  to  the  values  of  the  solution  on  a  specified  set  of  t  values. 
The  first  step  in  the  pointwise  solution  of  the  initial-value  problem,  dyjdt  =  f  {y,t)  with 
T  (^0 )  =  To  ’  is  to  approximate  the  value  of  y  at  h  =  tg  +  At ,  where  At  is  the  time  incre¬ 
ment.  This  new  value  of  y  is  generated  from  the  previous  value  of  y  plus  two  additional 
terms  that  are  proportional  to  the  derivative  evaluated  at  t  and  att  +  At ,  as  given  by 

Ti  =To+y/(t.0+y/(t.^+^0-  (14) 


After  has  been  generated  as  an  approximation  to  ^(ti)  the  procedure  can  be  repeated 

to  the  desired  length  of  t ,  which  is  time  in  this  simulation. 

B,  SYSTEM  DYNAMICS 

The  simulation  incorporates  a  fifth-order  binomial  series  representation  of  the  kill 
vehicle  guidance  system  components.  In  particular,  the  seeker  differentiates  the  line-of- 
sight  angle  and  is  represented  by  the  transfer  function 
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The  noise  filter  provides  a  smoothed  estimate  of  the  angular  line-of-sight  rate,  designated 
X ,  with  the  applicable  transfer  function  stated  as 
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Proportional  navigation  commands  are  generated  on  the  basis  of 


%  =  N'Vc.  (17) 
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Lastly,  the  relationship  between  aetual  and  commanded  acceleration  of  the  kill  vehicle  is 
dependent  on  the  flight  control  system,  which  is  represented  by  the  transfer  function 
given  by 
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The  flight  control  system  has  a  third-order  representation  to  accommodate  the  autopilot 
and  divert  thrusters,  which  comprise  the  divert  and  attitude  control  system. 

C.  SYSTEM  IMPLEMENTATION  IN  MATLAB 

A  schematic  representation  of  the  system  components,  including  the  variable 
names  used  in  the  MATLAB  code  is  presented  in  Figure  13.  All  transfer  functions  iden¬ 
tified  in  the  previous  section  are  implemented  in  code  using  the  variable  names  indicated 
in  Figure  13.  This  representation  of  the  individual  transfer  functions  permits  determina¬ 
tion  of  the  system  states  by  numerical  integration.  A  description  of  the  transformation  of 
Equation  9  into  the  components  of  loop  1  in  Figure  13  is  provided  in  Appendix  B. 


NT  Y  MISS  DISTANCE 


1  SEEKER 

2  NOISE  FILTER 

3-5  FLIGHT  CONTROL  SYSTEM 


Figure  13.  Simulation  of  Guidance  System  Components. 
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The  MATLAB  code  representing  the  kill  vehicle  guidance  system  is  attached  in 
Appendix  C.  The  simulation  facilitates  the  analysis  of  the  effects  of  various  acceleration 
saturations  and  how  they  impacts  miss  distance.  It  is  also  utilized  to  examine  how  an  ap¬ 
parent  target  displacement  during  the  terminal  maneuvers  affects  miss  distance.  In  both 
cases,  the  data  is  presented  in  a  normalized  form  to  ensure  transparent  summarization. 
Miss  distance  is  normalized  with  respect  to  displacement  at  initial  condition,  and  the  total 
flight  time  is  normalized  with  respect  to  guidance  system  time  constant.  The  system  time 
constant  is  a  critical  parameter  since  it  reflects  a  percentage  of  total  time  required  for  sys¬ 
tem  response  to  an  input  signal. 

1,  Acceleration  Saturation 

An  ideal  guidance  system  employing  proportional  navigation  will  always  reduce 
the  miss  distance  to  zero  if  the  kill  vehicle  has  sufficient  acceleration  capability.  The 
closed  form  solution  for  acceleration  required  by  a  kill  vehicle  with  a  zero  lag  guid¬ 
ance  system  engaging  an  accelerating  target  rij  is  given  by  [17] 


n 


c 


Hj 


N' 

N'-2 


f 

1 

V 


,N'-2 


‘■F  J 


(19) 


Even  for  this  idealized  situation,  a  nonzero  miss  distance  will  occur  if  the  kill  vehicle- 
target  maximum  acceleration  ratio  defined  by 


does  not  satisfy  the  inequality 
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In  Figure  14,  the  acceleration  ratio  is  plotted  for  the  three  common  values  of  N' , 

with  respect  to  normalized  time  t/fy  .  It  is  apparent  that  maximum  acceleration  is  re¬ 
quired  at  the  end  of  the  engagement,  and  this  maximum  requirement  decreases  as  N'  is 
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increased.  For  N'  =  3,  4,  and  5,  the  value  of  //  is  equivalent  to  3,  2,  and  1.67,  respee- 
tively. 


Figure  14.  Normalized  KV  Aceeleration  for  Zero  Lag  System. 


This  elosed  form  solution  is  not  applicable  to  the  more  realistic  guidance  system 
examined  herein;  therefore,  the  simulation  is  used  to  determine  the  value  of  ju  and, 
hence,  the  KV  aceeleration  capability  required  to  intercept  an  intermediate  range  ballistic 
missile,  which  has  an  average  aeceleration  of  3g,  a  bum  time  of  90  seconds,  and  a  range 
of  1000  kilometers  [1].  These  values  are  representative  of  the  North  Korean  No-Dong 
intermediate  range  ballistie  missile,  which  is  based  on  the  older  technology  Scud  missile. 
The  simulation  may  easily  be  altered  to  aecommodate  other  threat  missiles.  Engagement 
outcomes  are  plotted  for  a  range  of  values  from  2  to  infinity,  with  respect  to  miss 

distance  versus  total  flight  time  remaining  after  target  resolution.  The  target  lateral  ae¬ 
celeration  is  perceived  as  a  maneuver  by  the  KV  guidance  system;  therefore,  miss  dis- 
tanee  is  normalized  with  respect  to  KV-target  separation  distance  attributed  to  the  ma- 
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neuver.  The  simulation  is  run  separately  for  the  three  effective  navigation  ratios  of 
N'  =  3,  4,  and  5.  Graphs  are  generated  for  each  run  and  displayed  in  Figures  15,  16,  and 
17,  respectively.  In  each  case,  plots  for  infinite  acceleration  capability  are  shown  for 
comparison  purposes  only.  The  run  time  represents  the  total  simulated  engagement 
time,  which  commences  at  target  resolution  in  this  model,  and  is  designated  as  Total 
Flight  Time  in  the  following  graphs.  The  simulation  run  time  tp  is  not  to  be  construed  as 
the  total  flight  time  of  a  missile  from  boost  phase  to  terminal  phase.  The  simulation  also 
incorporates  a  time-loop,  which  progressively  increases  tp  to  examine  scenarios  of  in¬ 
creasing  duration  and  thereby  determine  the  effect  of  guidance  time  on  miss  distance. 


Figure  15.  Normalized  Miss  Distance  Due  To  Target  Acceleration  for  N'  =  3  . 

Figure  15  displays  the  simulation  results  for  a  kill  vehicle  guidance  system  utiliz¬ 
ing  N'  =  3.  It  is  clear  that  an  acceleration  ratio  of  2  or  less  will  fail  to  reduce  the  miss 
distance  to  zero.  Increasing  the  acceleration  ratio  to  3  also  fails  to  minimize  the  miss  dis- 


28 


tance  within  a  useful  time  frame.  The  minimum  acceleration  ratio,  which  effectively  re¬ 
duces  miss  distance  to  zero,  is  =  4.  This  indicates  that  a  kill  vehicle  with  accelera¬ 

tion  4  times  greater  than  the  target  will  have  the  capability  to  reduce  the  miss  distance  to 
zero  for  flight  times  greater  than  nine  times  the  guidance  system  time  constant.  In  the 
case  where  the  target  average  acceleration  is  3g,  the  KV  acceleration  required  would  be 
12g. 


Figure  16.  Normalized  Miss  Distance  Due  To  Target  Acceleration  for  A'  =  4 . 

Figure  16  shows  the  results  for  N'  =4,  where  the  miss  is  approximately  zero  for 
flight  times  that  are  ten  times  greater  than  the  guidance  system  time  constant,  and  the 
saturation  limit  is  at  least  3  times  the  target  acceleration.  Lastly,  Figure  17  shows  the  re¬ 
sults  for  N'  =  5,  where  the  miss  is  approximately  zero  for  flight  times  that  are  twelve 
times  greater  than  the  guidance  system  time  constant,  and  the  saturation  limit  is  at  least  2 
times  the  target  acceleration.  The  results  of  these  simulations  are  summarized  in  Table  2. 
It  is  apparent  that  increasing  the  effective  navigation  ratio  N'  results  in  decreasing  the 
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acceleration  requirement  for  the  kill  vehicle;  conversely,  the  total  flight  time  required  for 
intercept  increases. 


Figure  17.  Normalized  Miss  Distance  Due  To  Target  Acceleration  for  N'  =  5  . 


Effective  Navigation 

Ratio  (iV') 

Total  Flight  Time  After 

Target  Resolution  (/^) 

KV  Acceleration 

Requirement  ( ) 

3 

9t 

4nj^ 

4 

IOt 

Zhj 

5 

\2t 

liij 

Table  2.  KV  Acceleration  Requirements  for  Zero  Miss  Distance. 
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2,  Aimpoint  Displacement 

An  ideal  proportional  navigation  system  that  experiences  an  aimpoint  displace¬ 
ment  during  its  terminal  maneuvers  will  recover  and  ultimately  hit  the  sweet  spot,  assum¬ 
ing  adequate  acceleration  capability.  A  closed  form  solution  of  the  KV  acceleration  ca¬ 
pability  required  to  impact  a  target  after  a  step  displacement  of  the  aimpoint  is  given  by 
[17] 


N'{\-tltFY^^  y,c 


He  = 


(22) 


where  yj^  represents  initial  displacement  of  the  target  aimpoint.  Equation  22  states  that 
the  acceleration  requirement  is  greatest  at  the  time  when  aimpoint  displacement  occurs, 
and  that  is  proportional  to  the  amount  of  displacement  .  It  can  also  be  deduced 
that  miss  distance  is  zero  at  t  .  The  normalized  acceleration  provided  by  the  closed- 
form  solution  is  plotted  for  the  three  common  values  of  N'  in  Figure  18. 


Time/Total  Flight  Time  (t/tp) 


Figure  18.  Normalized  Acceleration  Due  To  Aimpoint  Shift  For  Ideal  System. 
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An  approximate  demarcation  point  can  be  seen  in  Figure  18  at  «  0.25  ;  up  to  that 

time,  increasing  the  effective  navigation  ratio  from  3  to  4  or  5  resulted  in  guidance  com¬ 
mands  necessitating  greater  acceleration  capability  to  execute,  and  subsequent  to  that 
time  the  acceleration  requirement  decreased  for  increasing  N' . 

A  closed-form  solution  is  not  available  for  determining  the  miss  distance  trig¬ 
gered  by  aimpoint  displacement  in  a  fifth-order  system,  as  used  in  this  simulation.  There¬ 
fore,  KV-target  engagement  outcomes  are  generated  on  the  basis  of  the  dimensionless 
parameter  {DP)  given  by 

0  5/2 

DP  =  — — ^252^  .  (23) 

yic 

The  parameter  DP  is  directly  proportional  to  the  maximum  acceleration  capabil¬ 
ity  of  the  kill  vehicle /2(,j^^^ .  Figures  19,  20,  and  21  display  normalized  miss  distance 

curves  attributed  to  an  aimpoint  displacement  equivalent  to  the  initial  condition  .  Sys¬ 
tems  with  effective  navigation  ratios  of  A'  =  3,  4,  and  5  are  considered  separately.  In 
each  case,  graphs  are  generated  with  DP  equivalent  to  0.1,  0.2,  0.5,  1,  and  infinity  for 
comparison.  As  in  the  previous  simulation,  each  increment  of  tp  represents  a  separate 
KV-target  engagement  time  frame  resulting  in  a  specific  miss  distance.  For  example,  at 
tp  =  0  there  would  be  zero  homing  time  remaining  after  the  aimpoint  displacement;  con¬ 
sequently,  the  kill  vehicle  would  miss  the  sweet  spot  by  an  amount  equivalent  to  the  total 
aimpoint  displacement  or  y  =  . 

Examination  of  Figures  19,  20,  and  21  reveals  that  acceleration  capability  is  in¬ 
versely  proportional  to  the  total  flight  time  or  number  of  system  time  constants  required 
to  minimize  the  miss  distance  to  zero.  In  other  words,  as  the  value  of  DP  is  increased, 
the  number  of  time  constants  required  is  decreased.  However,  kill  vehicles  with  large 
acceleration  capability  are  prone  to  increased  overshoot  of  the  target  and  may  reach  an 
unwanted  miss  distance  that  is  greater  than  the  original  aimpoint  displacement  as  evi¬ 
denced  in  Figure  21  for  DP  =  1  at  approximately  tp  =  4.5r .  Conversely,  there  is  a  prob¬ 
ability  for  kill  vehicles  with  less  acceleration  capability  to  incur  a  smaller  miss  distance 
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than  the  highly  energetic  kill  vehicles.  In  the  boost  phase  intercept  scenario,  zero  miss 
distance  is  ideally  required  to  ensure  complete  destruction  of  target  payload. 

Figure  19  shows  the  results  for  N'  =  3,  where  the  miss  distance  is  approximately 
zero  for  flight  times  that  are  eight  times  greater  than  the  guidance  system  time  constant. 
Figure  20  shows  the  results  for  N'  =  4,  where  the  miss  distance  is  approximately  zero  for 
flight  times  that  are  nine  times  greater  than  the  guidance  system  time  constant.  Lastly, 
Figure  21  shows  the  results  for  =  5,  where  the  miss  distance  is  approximately  zero  for 
flight  times  that  are  ten  times  greater  than  the  guidance  system  time  constant.  The  results 
of  this  simulation  are  summarized  in  Table  3.  Use  of  a  time  constant  less  than  or  equiva¬ 
lent  to  one-tenth  the  flight  time  remaining  after  aimpoint  displacement  will  ensure  the 
effectiveness  of  KV  guidance  system  response  for  all  three  values  of  N' . 


Figure  19.  Normalized  Miss  Distance  Due  to  Aimpoint  Displacement  for  N'  =  3  . 
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Miss  Distance/Aimpoint  Displacement  (yA^iq) 


Figure  20.  Normalized  Miss  Distance  Due  to  Aimpoint  Displacement  for  A'  =  4 . 


Figure  2 1 .  Normalized  Miss  Distance  Due  to  Aimpoint  Displacement  for  N'  =  5  . 
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Effective  Navigation  Ratio  (iV' ) 

Total  Flight  Time  (/^)  Remain¬ 
ing  After  Aimpoint  Displacement 

3 

8r 

4 

9r 

5 

lOr 

Table  3.  Aimpoint  Displacement  Results. 

The  graphs  displayed  in  Figures  19,  20,  and  21  may  be  used  to  determine  the 
amount  of  miss  distance  incurred  for  any  displacement  of  the  aimpoint  when  the  maxi¬ 
mum  acceleration  capability  and  guidance  system  time  constant  of  a  kill  vehicle  are 
specified.  In  our  case,  the  KV  acceleration  requirements  generated  in  the  previous  sec¬ 
tion  are  used  in  conjunction  with  the  prerequisite  for  zero  miss  distance  to  determine  a 
suitable  time  constant  for  the  KV  guidance  system.  Equation  23  is  employed  to  calculate 
the  time  constants  for  DP  =  0.1  and  DP  =  1.0  curves,  listed  in  Table  4.  The  original 
aimpoint  displacement  is  assumed  to  be  equivalent  to  60  feet,  based  on  the  length  of  the 
No-Dong  missile  plus  an  offset  for  the  location  of  the  power  centroid  of  the  plume. 


DP 

Aiij. 

IHj. 

0.1 

0.18 

0.20 

0.25 

1.0 

0.56 

0.64 

0.79 

Table  4.  KV  Guidance  System  Time  Constants  (Seconds). 

To  minimize  the  probability  of  overshoot,  the  DP  =  0.1  curve  is  most  applicable; 
secondly,  to  reduce  the  intercept  time,  it  is  necessary  to  use  the  guidance  system  with 
N'  =  3  ,  which  requires  a  maximum  acceleration  capability  of  =  Arij .  On  the  basis  of 
these  considerations,  it  is  recommended  that  a  time  constant  of  0.18  seconds  be  utilized. 
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Use  of  r  =  0.18  seconds  in  conjunction  with  tp  =  lOr  seconds  enables  the  kill  vehicle  to 
reduce  the  miss  distance  to  zero  for  an  aimpoint  displacement  occurring  up  to  1.8  seconds 
prior  to  termination  of  the  engagement. 

The  simulation  facilitated  the  analysis  of  acceleration  saturation  and  aimpoint  dis¬ 
placement  in  a  proportional  navigation  guidance  system.  Minimum  requirements  for  the 
kill  vehicle  acceleration  capability  and  guidance  system  time  constant  were  determined  to 
ensure  a  negligible  miss  distance  at  termination  of  engagement.  The  results  are  summa¬ 
rized  in  the  concluding  chapter. 
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V.  CONCLUSION 


Boost  phase  interception  of  ballistic  missiles  is  constrained  by  the  dual  require¬ 
ments  of  minimizing  interception  time  to  preempt  deployment  of  submunitions  or  multi¬ 
ple  reentry  vehicles  and  reducing  miss  distance  to  within  a  few  centimeters  of  the  war¬ 
head  sweet  spot  to  ensure  complete  destruction  of  target  payload. 

The  simulation  implements  proportional  navigation  commands,  in  conjunction 
with  a  fifth-order  binomial  series  representation  of  the  system  dynamics,  to  assess  the 
impact  of  acceleration  saturation  on  miss  distance.  The  simulation  was  also  employed  to 
assess  the  influence  of  aimpoint  displacement  on  miss  distance  during  the  terminal  ma¬ 
neuvers  of  the  engagement.  The  target  characteristics  utilized  were  that  of  a  No-Dong 
intermediate  range  ballistic  missile. 

The  analysis  of  acceleration  saturation  revealed  that  increasing  the  effective  navi¬ 
gation  ratio  N'  results  in  decreasing  the  acceleration  requirement  for  the  kill  vehicle  to 
attain  zero  miss  distance.  Conversely,  decreasing  N'  reduces  the  total  flight  time  re¬ 
quired  for  interception.  Since  time  is  the  critical  factor  for  boost  phase  interception  of 
ballistic  missiles,  it  is  recommended  that  a  guidance  system  with  N'  =  3  be  implemented 
for  the  kill  vehicle  guidance  system.  This  results  in  a  12g  lateral  acceleration  require¬ 
ment  for  the  kill  vehicle  to  completely  destroy  a  No-Dong  missile. 

The  analysis  of  aimpoint  displacement  confirmed  that  the  acceleration  require¬ 
ment  is  inversely  proportional  to  the  flight  time  necessary  to  minimize  the  miss  distance 
to  zero.  Unfortunately,  kill  vehicles  with  large  acceleration  capability  are  prone  to  in¬ 
creased  overshoot  of  the  target  and  may  reach  an  unwanted  miss  distance  that  is  greater 
than  the  original  aimpoint  displacement.  The  optimized  time  constant  for  the  simulated 
guidance  system  was  determined  to  be  0.18  seconds  or  one-tenth  the  flight  time  remain¬ 
ing  after  aimpoint  displacement. 

Recommendations  for  future  research  include  determination  of  the  following  pa- 
rametrics: 

•  Distance  from  which  the  onboard  sensor’s  optics  can  identify  the  target 
payload; 
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•  Fuel  capacity  required  onboard  the  kill  vehicle  to  sustain  the  divert  rocket 
thrusters;  and 

•  Maximum  velocity  required  to  perform  divert  maneuvers. 
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APPENDIX  A.  ALPHA  KV-MODEL 


A  listing  of  the  MATLAB  commands  used  to  implement  the  ALPHA-KV  model 
is  provided  in  this  appendix. 

% - 

LT=0.9;  %Target  Payload  Length  (m) 

DT=0.6;  %Target  Payload  Diameter  (m) 

DS=0.08;  %Target  Submunition  Diameter  (m) 

rhot=1.5;  %Density  of  Target  (g/cm^3) 

rhokv=1.0;  %Density  of  Kill  Vehicle  (g/cm^3) 

BHN=25;  %Target  Brinell  Hardness  Number 

co=9.5e4;  %Sound  Speed  of  Target 

NT=3;  %Number  of  Sub  munition  Tiers 

epsilon=0.03;  %Required  Overlap  Distance  of  Crater  into  Last  Tier 

alpha=3.54e-10; 

Lc=(LT*((NT-l)/NT))+epsilon;  %Required  Crater  Depth 
%Minimum  Radius  to  Overlap  All  Submunitions 
Ro=sqrt((Lc-0.5*LT)^2+(0.5*DT-DS)^2); 

Rm=0.0;0. 1:0.3; 
vc=pi*Lc*(Ro+Rm)^2; 

V=L5:0.05:6; 
for  rm=l  :length(Rm) 
for  v=l:length(V) 

Lkv(v)=Lc*sqrt(rhot/rhokv);  %KV  Length  (m) 

m(rm,  v)=(2  *pi  *  (Lc  *  Ro^2+Lc  *  Rm(rm)  *  (2  *  Ro+Rm(rm))))  *  0 . 1 /(alpha*  ( V(v)  *  1 00 

0r2); 

Dkv(rm,v)=100*sqrt(4*m(rm,v)*0.001/(pi*rhokv*Lkv(v))); 

E(rm)=pi  *Lc  *  (Ro^2+Rm(rm)  *  (2  *  Ro+Rm(rm)))  *  0 . 1  /alpha; 

end 

end 

plot(V,m') 

grid 

legend('R_m  =  0','R_m  =  10cm','R_m  =  20cm', 'R_m  =  30cm') 
title(['Optimal  Kill  Vehicle  Mass  vs  Relative  Velocity']) 
ylabel('Minimum  KV  Mass  (kg)') 
xlabel('Impact  Velocity  (km/s)') 

% 

figure(2) 

plot(V,Dkv',V,100*Lkv,'k-') 

legend('R_m  =  0','R_m  =  10cm','R_m  =  20cm', 'R_m  =  30cm','KV  Length') 
title(['Optimal  Kill  Vehicle  Diameter  vs  Relative  Velocity']) 
ylabel('KV  Diameter  (cm)') 
xlabel('Impact  Velocity  (km/s)') 
grid 
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% 

figure(3) 

plot(Rm*  1 00, E/ 1  e6) 
grid 

%title(  ['Minimum  Kinetic  Energy  for  Optimized  KV']) 
ylabel('Required  Kinetic  Energy  (MJ)') 
xlabel('Aimpoint  Offset  (cm)') 
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APPENDIX  B.  TRANSFER  FUNCTION  DERIVATION 


The  individual  transfer  funetions  shown  in  Figure  1 1  are  transformed  into  the  loop 
eomponents  of  Figure  13,  to  faeilitate  determination  of  the  system  states  via  numerieal 
integration.  The  seeker  bloek  is  examined  here  and  is  representative  of  all  the  other 
transfer  funetion  bloeks.  For  eonvenienee  Figure  22  shows  the  equivalenee  of  eompo- 
nents. 


SEEKER  REPRESENTATION  IN  FIGURE  11 


SEEKER  REPRESENTATION  IN  FIGURE  13 


JL 


Figure  22.  Transfer  Funetion  Representation. 


The  variable  definitions  for  the  above  eomponents  are  as  follows: 

•  A  in  Figure  1 1  is  defined  as  LAMDA  in  Figure  13  and  in  the  simulation; 

•  A  in  Figure  1 1  is  defined  as  LAMDA  RATE  in  Figure  13  and  in  the 

simulation; 

•  5  in  Figures  1 1  and  13  is  the  Laplaee  operator; 

•  SD  in  Figure  13  is  an  internal  state  of  the  seeker;  and 

•  r  in  Figures  1 1  and  13  is  defined  as  TAU  in  the  simulation. 

The  MATLAB  eode  representing  the  seeker  output  is  determined  from  the  right- 
hand  diagram  in  Figure  22,  as 

LAMDA  RATE  =  ^^^(LAMDA  -  SD) .  (24) 
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To  prove  Equation  24  is  directly  related  to  the  seeker  transfer  function  shown  in  the  left- 
hand  diagram  in  Figure  22,  note  that 


SD  _  1 
LAMDA  RATE  ^  5 


(25) 


and  substitute  Equation  25  into  24,  as  follows 


EAMDA  RATE  = 


5 

TAU 


( 

LAMDA 

V 


LAMDARATE^ 

) 


(26) 


LAMDA  RATE*TAU  =  5*LAMDA-5^ 


LAMDA  RATE 


LAMDA  RATE 


^  5^ 

TAU  +  - 

s) 


=  5*LAMDA 


LAMDA  RATE 


LAMDA 

V  s) 


(27) 


(28) 


(29) 


LAMDARATE 

LAMDA 


5 

"5*TAU  +  5^ 
V  J 


(30) 


LAMDA  RATE  _  5^ 

LAMDA  ~5*TAU  +  5 


(31) 


LAMDARATE 

LAMDA 


s 

(,  5*TAU^ 

I  5  ] 


(32) 
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Equation  32  may  be  re-written  with  the  symbols  utilized  in  the  left-hand  diagram  of  Fig¬ 
ure  22  as 


A 


s 


^  st'^ 
1+— 


which  is  the  original  transfer  function. 


(33) 
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APPENDIX  C.  GUIDANCE  SYSTEM  SIMULATION 


A  listing  of  the  MATLAB  commands  used  to  implement  the  fifth-order  propor¬ 
tional  navigation  guidance  simulation  is  provided  in  this  appendix. 


%Variable  Definitions — 

%NT 

%YDD 

%YD 

%Y 

%LAMDA 

%LAMDA_RATE 

%SD 

%ESTDD 

%EST_EAMDA_RAT 

%NC 

%X4D 

%X4 

%X5D 

%X5 

%NLD 

%NL 

%Begm  Simulation - 

VC=32808. 

G=32.2 

NT=3*G 

DISPEACE=1. 

TAU=1. 

NPRIME=4. 

NCEIM=10el0 

n=0 

for  TP=0. 1:0. 1:12 

Y=DISPEACE 

YD=0. 

NE=0. 

SD=0. 

EST_EAMDA_RATE=0. 

X5=0. 

T=0. 

DEETA=.01 

while  T<=(TP-le-5) 

YOED=Y 

YDOED=YD 

NEOED=NL 

SDOED=SD 


Target  Acceleration 

Relative  Acceleration 

Relative  Velocity 

Miss  Distance 

Geometric  EOS 

Geometric  EOS  Rate 

Seeker  Internal  State 

Eilter  Internal  State 

Eilter  Estimate  Geometric  EOS  Rate 

Commanded  Acceleration 

ECS  Autopilot  Internal  State 

ECS  Autopilot  Output 

ECS  Divert  Internal  State 

ECS  Divert  Output 

ECS  Internal  State 

Achieved  Acceleration 


YoClosing  Velocity 

%Acceleration  due  to  gravity 

%Target  Eateral  Acceleration 

%Aimpomt  Displacement 

%Kill  Vehicle  Guidance  System  Time  Constant 

%Effective  Navigation  Ratio 

%Acceleration  Saturation  Infinite  Limit  for  KV 

%Pass  Number 

%Outer  Loop  for  Generating  Total  Elight  Time 

%Initial  Condition  Displacement 

%Enter  Iteration  Process  Velocity 

%Enter  Iteration  Process  Acceleration 

%Enter  Iteration  Process  LOS  Angle 

%Enter  Iteration  Process 

%Enter  Iteration  Process 

%Begm  Iteration  Process  at  Time  =  0 

%Increment  Length  for  R-K  Numerical  Integration 

%Run  Iteration  Until  TGO  Less  Than/Equal  le-5  s 

%Miss  Distance 

%Miss  Distance  Derivative 

%Acceleration  Achieved  by  Kill  Vehicle 

%Output  from  Seeker  Integration 
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EST_LAMDA_RATE_OLD=EST_LAMDA_RATE  %Estimate  EOS  Rate 

X40ED=X4  %State  4  within  Elight  Control  System 

X50ED=X5  %State  5  within  Elight  Control  System 

STEP=1 

EEAG=0 

while  STEP<=1 

if  PEAG==1  %Conduet  Iteration  at  Time  T 

STEP=2 

Y=Y+DEETA*YD  %Displacement  at  Time  T 

YD=YD+DEETA*YDD  %Derivative  of  Position  at  Time  T 

NE=NL+DEETA*NLD  %Achieved  Acceleration  at  Time  T 

EST_EAMDA_RATE=EST_LAMDA_RATE+DEETA*ESTDD 

SD=SD+DEETA*LAMDA_RATE 

X4=X4+DEETA*X4D 

X5=X5+DEETA*X5D 

T=T+DEETA  %Increment  Time 

end 

%Conduct  Iteration  at  Time  T+DEETA - 

TGO=TE-T+O.OOOOI  %Time  To  Go 

EAMDA=Y/(VC*TGO)  %Measured  Geometric  EOS  Angle 

EAMDA_RATE=5.0*(EAMDA-SD)/TAU  %Seeker  Output  EOS  Rate 

%Eilter  Internal  Estimated  Geometric  EOS  Rate 
ESTDD=5.0*(LAMDA_RATE-EST_EAMDA_RATE)/TAU 

%Commanded  Acceleration  from  PN  Eaw 
NC=NPRIME*VC*EST_LAMDA_RATE 

ifNC>NCEIM 

NC=NCEIM  %Acceleration  Eimit  of  KV 

end 

ifNC<-NCEIM 

NC=-NCEIM  %Acceleration  Eimit  of  KV 

end 

X4D=5.0*(NC-X4)/TAU  %Ehght  Control  System 

X5D=5.0*(X4-X5)/TAU  %Elight  Control  System 

NED=5.0*(X5-NE)/TAU  %  Achieved  Acceleration  of  KV 

YDD=NT-NE  %Relative  Acceleration 

EEAG=1 
end 

EEAG=0 

Y=0.5*(YOED+Y+DEETA*YD)  %Runge-Kutta  Integration 

YD=0.5*(YDOED+YD+DEETA*YDD)  %Runge-Kutta  Integration 

NE=0.5*(NEOLD+NL+DELTA*NLD)  %Runge-Kutta  Integration 

SD=0.5*(SDOED+SD+DELTA*EAMDA_RATE)  %Runge-Kutta  Integration 
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EST_LAMDA_RATE=0.5*(EST_LAMDA_RATE_OLD+EST_LAMDA_RATE 
+DEETA*ESTDD)  %Runge-Kutta  Integration 

X4=0.5*(X4OLD+X4+DEETA*X4D);  %Runge-Kutta  Integration 

X5=0.5*(X5OED+X5+DEETA*X5D);  %Runge-Kutta  Integration 

end 

n=n+I  %Inerement  Simulation  Pass  Number 

Array  TE(n)=TE 

Array  Y(n)=Y 

end 

figure 

plot(  ArrayTE,  Array  Y,'r—') 
grid 

xlabel('Elight  Time  (See)') 

ylabel('Miss  (Et)') 

ele 

output=[ArrayTE',Array  Y'] ; 
save  datfil.txt  output  -ascii 
disp  'simulation  finished' 
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